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ABSTRACT 

We examine diffusive shock acceleration (DSA) of the pre-exisiting as well 
as freshly injected populations of nonthermal, cosmic-ray (CR) particles at weak 
cosmological shocks. Assuming simple models for thermal leakage injection and 
Alfvenic drift, we derive analytic, time-dependent solutions for the two popu- 
lations of CRs accelerated in the test-particle regime. We then compare them 
with the results from kinetic DSA simulations for shock waves that are expected 
to form in intracluster media and cluster outskirts in the course of large-scale 
structure formation. We show that the test-particle solutions provide a good 
approximation for the pressure and spectrum of CRs accelerated at these weak 
shocks. Since the injection is extremely inefficient at weak shocks, the pre-existing 
CR population dominates over the injected population. If the pressure due to 
pre-existing CR protons is about 5 % of the gas thermal pressure in the upstream 
flow, the downstream CR pressure can absorb typically a few to 10 % of the shock 
ram pressure at shocks with the Mach number M < 3. Yet, the re-acceleration of 
CR electrons can result in a substantial synchrotron emission behind the shock. 
The enhancement in synchrotron radiation across the shock is estimated to be 
about a few to several for M ~ 1.5 and 10^ — 10^ for M ~ 3, depending on 
the detail model parameters. The implication of our flndings for observed bright 
radio relics is discussed. 

Subject headings: acceleration of particles — cosmic rays — galaxies: clusters: 
general — shock waves 
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INTRODUCTION 



Cosmological shock waves result from supersonic flow motions in duced by hierarchi- 



cal clustering du ring the large-scale structure formation in the Universe (IMiniati et al.ll2000 



Ryu et al.ll2003[ ). According to studies based on cosmological hydrodynamic simulations, the 
shocks formed by merger of subclumps, infall of matter and internal flow motion in intraclus- 



ter media 



ICMs ) and cluster outskirts are relatively weak with Mach number M < a few 



(iRvu et al. 



20081: 



12003 



Vazza et al. 



frommer et al.l l2006t iKang et al.l 120071 : ISkillman et al. 



20091). Indeed, ob servations of X-ray s hocks (e.g., 



20101: 



van Weeren et al. 



2005 : lMarkevitch fc Vikhlinir]j|2007n and radio relics (e.g.. lBagchi et al.ll2006uFinoguenov et al 



2008 



Hoeft et al. 



Markevitch et al. 



2002 



20101 ) indicate that the estimated Mach number of observed shocks 



in cluster environments is consistent with such theoretical predictions. 

Suprathermal particles are known to be produced as an inevitable consequence of the 
formation of collisionless shocks in tenuous plasmas and they can be further accelerated to 
become cosmic rays (CRs) through interactions with resonantly scattering Alfve n waves in 
the converging flow across a shock ( lBelllll978l : lDrurylll983l : iMalkov fc Druryll200ll ). Detailed 
nonlinear treatments of diffusive shock acceleration (DSA) have predicted that at strong 
shocks a significant fraction of the shock kinetic energy is transfer red to CRs, inducing 



highly nonlinear back-reactions from CRs to the underlying flow (e.g., lAmato fc Blasil 12006 



Vladimirov et al.l l2006l : iKang fc JonesI 120071 ) . Multi-band observations of nonthermal radio 



Parizot et al. 


2006; 


Reynolds 


2008; 


Abdo et al. 


2010) 



The presence of nonthermal particles, especially electrons, in clusters of galaxies, has 
been inferred from observations of synchrotron em ission from radio halos and relics (see. 



e.g., iCarilli fc Taylorl l2002t iGovoni fc Ferettil l2004j . for review). Since the matter in ICMs 
and cluster outskirts should have gone first through accretion shocks of high Mach number 
around nonlinear structures and th en through weaker shocks due to mergers and flow motion 
( jRyu et al.ll2003l : iKang et al.ll2007l ). DSA should be responsible for at least a part of the CR 
production. Nonthermal particles can be also produce d via turbulent acceleration (see, e.g., 
Cassano fc Brunettil 120051 : iBrunetti fc LazarianI 12007 ). Recent Fermi observations of 7-ray 
emission from galaxy clusters, however, lim it that the pressure due to CR proton s cannot 
exceed ~ 10 % of the gas thermal pressure ( lAbdo et al.ll2010t iDonnert et al.ll2010l ). 



At weak shocks with M < a few, DSA is known to be rather inefficient and the CR 
pressure remains dynamically insignifi cant , partly beca use the injection from thermal to 
nonthermal particles is inefficient (e.g., iKang et al.l 120021 ). In such test-particle regime, the 



downstream CR spectrum takes the power-law form of f^iv) oc p where the spectral 



slope, g, depends on the velocity jump across the shock ( lDrurylll983l ). Recently, iKang fc Ryu 
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( 120101 ) suggested analytic, time-dependent solutions for the test-particle CR spectrum, using 
results from DSA simulations in which particles are injected via thermal leakage process and 
accelerated to ever increasing maximum momentum, Pmax(^)- They found that at weak shocks 
expected to form in ICMs and cluster outskirts, indeed, much less than ~ 10~^ of particles 
are injected into CRs and much less than ~ 1% of the shock ram pressure is converted into 
the downstream pressure of CR protons, so the particle acceleration is virtually negligible. 

However, the rec ent discovery of very bright radio relics associated with weak shocks 



of M < a few (e.g., iBagchi et al.l l2006l : iFinoguenov et al.l I2OIOI : Ivan Weeren et al.l I2OIOI ) 



suggests that, contrary to the expectation, DSA should operate at weak shocks in cluster 
environments. One way to explain this is to presume that the relics form in media with 
pre-existing CRs which were produced by DSA at previous shocks and/or by turbulent 
acceleration. The existence of pre-exiting CRs alleviates the problem of inefficient injection 
at weak shocks. 

In this paper, we examine the DSA at weak cosmological shocks in the presence of 
pre-existing CRs. First, the properties of weak shocks in ICMs and cluster outskirts are 
briefly reviewed in Section 2. Analytic, time-dependent solutions for the acceleration of the 
pre-existing and freshly injected populations of CRs in the test-particle regime is described 
in Section 3, while the numerical solutions from kinetic DSA simulations are presented 
in Section 4. The synchrotron radiation from CR electrons accelerated at these shocks is 
discussed in Section 5. Finally, a brief summary is given in Section 6. 



SHOCK WAVES IN ICMS AND CLUSTER OUTSKIRTS 



Shock waves in the large-scale structure of the universe have been studied in details using 
various hydrodyna mic simulations for the cold dark matter cosmology wi t h cosmological 



constant (ACDM) fiRvu et al. 



2008: Hoeft et al 



200^ 



12003 



Pfrommer et al.ll2006l : iKang et al 



2007: Skillman et al. 



sT vazza et al. 2009 ). It was found that shocks with Mach number 



typically up to M ~ 10^ and speed up to ~ a few xlOOO km s~^ at the present universe 
{z = 0). In ICMs and cluster outskirts, however, shocks are expected to have lower Mach 
number, because they form in the hot gas of kT > keV. 

To examine the characteristics of shocks in ICMs and cluster outskirts, we analyze the 
shocks with the preshock gas temperature of Ti > 10'' K. The c osmic web is fille d with 
ionized plasmas, the intergalactic medium (jCen &: Ostrikerl Il999l : iKang et al.ll2005l ). The 
hot gas with T > 10'' K is found mostly in ICMs and cluster outskirts, and the Warm Hot 
Intergalactic Medium (WHIM) with 10^ K < T < 10'' K is distributed mostly in filaments. 
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The diffuse gas with T < 10^ K resides mainly in sheetlike structures and voids. The 
shocks were found in a simulation of the WMAPl-normalized ACDM cosmology employed 
the following parameters: = 0.048, fi„ = 0.31, = 0.69, h = //o/(100 km/s/Mpc) 
= 0.69, (Tg = 0.89, and n = 0.97. The simul ation was performed using a PM/Eulerian 
hydrodynamic cosmology code (IRyu et al.lll993[ ) . Detailed descri p tions for numerical set-up 
and input physical ingredients can be fo und in ICen fc Ostrikerl (120061 ). The procedure to 
identify shocks was described in details in Ryu et al. ( 20031 ). 

Figure 1 shows the surface area of shocks with Ti > 10'' K per Mach number interval 
in the entire simulation volume, normalized by the volume. Here, S is given in units of 
{h~^Mpc)~^ . The quantity 5* provides a measure of shock frequency or the inverse of the 
mean comoving distance between shock surfaces. To avoid confusion from complex flow 
patterns and shock surface topologies associated with very weak shocks, only those portions 
of shock surfaces with M > 1.5 are shown. We also calculated the incident shock kinetic 
energy flux, = (l/2)pi'u^, where pi is the preshock gas density, and then the kinetic energy 
flux through shock surfaces per Mach number interval, normalized by the simulation volume, 
dF^{M) / dM . Figure 1 shows dF^{M)/dM, too. As expected, the Mach number of the shocks 
formed in ICMs and cluster outskirts is small, typically M < 3. The frequency increases 
to weakest possible shocks with M ~ 1. The kinetic energy flux through shock surfaces is 
larger for weaker shocks; that is, weaker shocks process more shock energy, conflrming the 
energetic dominance of weak shocks in cluster environments. 



ANALYTIC TEST-PARTICLE SPECTRUM 



In the kinetic DSA approach, the following diffusion-convection equation for the pitch- 
angle-averaged distribution function of CRs, f{x,p, t), is solved along with suitably modifled 
gasdynamic equations: 



dt 



+ [u + u^ 



dx 



p d{u + 
3 dx 



dp dx 



k{x,p) 



df 

dx 



(1) 



where n{x,p) is the spatial diffusion coefficier it and u,,, is th e drift speed of the local Alfvenic 
wave turbulence with respect to the plasma (jSkillinglll975l ). The scattering by Alfven waves 
tends to isotropize the CR distribution in the wave frame, which may drift upstream at the 
Alfven speed, va, with respect to the bulk plasma. So the wave speed is set to be Uw = —va 
upstream of shock, while = Q downstream. 

In the test-particle regime where the feedback due to the CR pressure is negligible, the 
downstream CR distribution can be described with a power-law spectrum, f2{p) oc and 
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the slope is given by 



3{ui - va) 

Ui-Va- U2 



3(t(1 



1 



(2) 



where Ui and U2 are the upstream and downstream flow speeds, respectively, in the shock 
rest frame, a = Ui/u2 = P2/P1 is the shock compr ession ratio, and Ma = U}Iva is the 
upstream Alfven Mach number with va = Bi/ yJAirpi (jDrurylll983l : iKang fc Ryijl2010l ). The 



test-particle power-law slope q can be calculated as a function of shock Mach number M 
with a = [(7g + l)M2]/[(7g - 1)M^ + 2], which becomes AM^/{M'^ + 3) for a gas adiabatic 
index 7^, = 5/3, and Ma = M/5. Here, 5 = va/cs is the Alfven speed parameter, where Cg is 
the upstream sound speed. The maximum momentum of CR protons achieved by the shock 
age of t can be estimated as 



Pmax(^) 



nipC 



:i-M7^ 



)(a-l 



M 



3a(2 - MJ^) 



K* 



(3) 



(4) 



where Ug = ui is the shock speed (jPruryl Il983l : iKang fc Ryull2010l ). Here, a Bohm-type 
diffusion coefficient, 

' P WPo 
jipc) V P 

is adopted, where k* = mpC^/{3eBo) = 3.13 x 10^^(i?o/l pG)~^cm^s~^, Bq and po are 
magnetic field strength and the gas density far upstream. In CR-modified shocks where CRs 
are dynamically non-negligible, in general, the upstream flow is decelerated in the precursor 
before it enters the gas subshock. So we use the subscripts "0", "1", and "2" to denote 
the conditions far upstream, immediate upstream and downstream of shock, respectively. 
Of course, in the test-particle limit, the distinction between far and immediate upstream 
quantities disappears, e.g., Po = Pi- 



In the limit of large M {a ? 
protons can be approximated by 



4) and large Ma {S ~ 0), the maximum energy of CR 



uH 



max,p 



-nipC 



10 



10 



GeVf- 



t 



Bo 



(5) 



VlO^kms-iy yiO^jTS J \1 fiG ^ 

The CR proton spectrum limited by the shock age is expected to have a cutoff at around 
~ Pmax(^) (see Section 3.3 for further discussion). 



3.1. Pre-existing Population 



As noted in Introduction, it seems natural to assume that ICMs and cluster outskirts 
contain pre-existing CRs. But their nature is not well constrained, except that < O.lPg, 
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i. e. , the pressure of CR protons is l ess that ~ 10 % of the gas thermal pressure (e.g., 
Abdo et al.ll2010l : iDonnert et al.ll2010l ). With pre-existing CRs of spectrum fo{p) upstream 
of shock, the steady-state, test-particle solution of Equation ([T]) for the downstream CR 
distribution can be written as 



Pinj 



VPinj 



(6) 



where q is the test-particle power-law slope given in Equation ([2]) (lDrurylll983l ). Here, pinj is 
the lowest momentum boundary above which particles can cross the shock, i.e., the injection 
momentum (see the next subsection). By this definition of Pinj, the CR distribution function, 
/o = and /2 = for p < pinj. The first term in the right-hand-side of Equation represents 
the re-accelerated population of pre-existing CRs, while the second term represents the 
population of CRs freshly injected at the shock and will be discussed in the next subsection. 

We adopt a power-law form, fo{p) = /pre ■ {p/Pm})~^, with the slope s = 4 — 5, as the 
model spectrum for pre-existing CR protons. If pre-existing CRs were generated at previous 
shocks, the slope of s = 4 — 5 is achieved for M > \/5 with 6 = (see Equation ([2])). On 
the other hand, if they are mainly the outco me of turbulent acceleration, the slope should 



be close to s ~ 4 (see, e.g., IChandranI l2005l ) . Then, the spectrum of re-accelerated CRs is 
obtained by direct integration: 



fr%p) 

liqj^s, for p >]9inj, 



[q/{q - s)] [1 



frip) 



\q 



(pMnj)-'^+1/o(p), 

gln(p/pinj)/o(p). 



-f 

I J pre 1 
S\ VPinj 



if g 7^ s 
if g = s. 



(7) 



where r = min(g,s). That is, if the spectral slope of pre-existing CRs is softer than the 
test-particle slope {s > q), the re-accelerated CR spectrum gets flattened to p~'^ by DSA; in 
the opposite case {s < q), the re-accelerated CR spectrum is simply amphfied by the factor 
of q/{q — s) and retains the same slope as the slope of pre-existing CRs. 

Figure 2 shows the re-accelerated CR distribution given in Equation ([7]) for a M = 3 
shock in the presence of the pre-existing power-law CR spectrum with the slope s = 4 and 4.5 
(right panel) and s = 5 (left panel). The Alfvenic drift is ignored (5 = 0), so the test-particle 
slope is g = 4.5. Here, we adopted the following parameters: the upstream gas temperature 
To = 10^ K and the injection parameter = 0.25, resulting in pinj = 8.0 x 10"'^ 
next subsection for details of our injection model). 



rUpC (see the 



The figure illustrates that for p ^ pinj, the CR amplification factor, f2{p)/ fo{p), ap- 
proaches a constant, q/{q — s) = 9, in the case of s = 4, increases as ln(p/pinj) in the case of 
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q = s = 4.5, and scales as (p/pinj)*"'' in the case of s = 5. So, for instance, the factor becomes 
/2//0 = 32 and 310 at p/nipC = 10 for s = 4.5 and 5, respectively. We point that these values 
of the CR amplification factor are substantially larger than those expected for the adiabatic 
compression across the shock. With pre-existing CRs of /o oc p~^, the amplification factor 
due to the adiabatic compression is given by 

/|*//o = cr'^' (9) 

in the test-particle regime. So the adiabatic amplification factor is /I'^^/Zo = 4.3, 5.2, and 
6.2 and for s = 4, 4.5 and 5, respectively, at a Mach 3 shock. Note that the adiabatic 
compression does not change the slope of the CR spectrum. 

The left panel of Figure 2 also shows the time evolution of the CR distribution at the 
shock location, fs{p, t), from a DSA simulation for the same set of parameters (see Section 4 
for details of DSA simulations). The CR injection was turned off for this particular simulation 
in order to compare the analytic and numerical solutions only for pre-existing CRs. This 
demonstrates that the time-dependent solution asymptotes to the steady-state solution in 
Equation ([7]). 



3.2. Injected Population 



Because complex plasma interactions among CRs, waves, and the underlying gas flow 
are not fully understood yet, it is not possible t o make a precise quanti tative prediction for 



Malkov fc Drury 



200 ll ). Here, we adopt a 



the injection process from first principles (e.g., 
phenomenological injection model that can emulate the thermal leakage process, through 
which particles above a certain injection momentum P i pj cro ss the shock and get injected 
to the CR population (IKang et al.l l2002l : iKang fc Ryull2010l ). Then, the CR distribution 
function at Pinj is anchored to the downstream Maxwellian distribution as 



/inj = /(Pinj) = ^ Pth exp (-Qfnj) . (10) 

where n2 is the downstream proton number density. Here, Pinj and Qinj are defined as 

1.07' 



Qi.j(M) ^ ^ ^ 1.17^ 
Pth Pth 



1 + 



M 

y 



0.1 



11 



where pth = A/2mJ^AJ^T2 is the thermal peak momentum of the downstream gas with tem- 
perature T2 and ks is the Boltzmann constant. We note that the functional form of Qinj was 
adopted to represent an "effective" injection momentum, since particles in the suprathermal 
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tail c an cross the shock with a smoothly-varying probability distribution (see iKang et al. 
20021 ). One free parameter that controls the leakage process is the injection parameter, 
eB = Bq/Bj_, which is the ratio of the general magnetic field along the shock normal, Bq, 
to the amplitude of the downstream, magnetohydrodynamic (MHD) wave turbulence, i?^. 
Although plasma hyb rid simulations and theories both suggested that 0.25 ^ < 0.35 
(IMalkov fc Volklll998l ). the physical range of this parameter remains to be rather uncertain 
due to lack of full understanding of relevant plasma interactions. 

The second term in Equation is fixed by q, Pinj, and /inj. The fraction of particles 
injected into the CR population can be estimated analytically as well: 



n2 



Qfnj exp (-Qfnj) 



(12) 



which is fixed only by Qmj and q. The inj ection fraction depe nds strongly on (through 
Qinj) for weak shocks with M < 5 (see also lKang fc Ryull2010l ). For example, it varies from 
5 X 10-^ to 10-3 for eb = 0.25 - 0.3 for shocks with M = 3. 



3.3. Cosmic-Ray Spectrum for Weak Shocks 



Kang fc Ryul ( l2010l ) demonstrated that the time- dependent, test-particle solutions of the 
downstream CR distribution can be represented by the steady-state, test-particle solutions 
with an exponential cutoff ( ICaprioli et al.l 120091 ). if the cutoff momentum is set as p* ~ 
1.2 Pmax(^) with Pmax(^) in Equatiou ([3]). Here, we suggest that the same cutoff would be 
applied to the spectrum of re-accelerated CRs. Then, the CR distribution at the shock 
location, Xg, originated from both the pre-existing and freshly injected populations can be 
approximated by 



fs{,P,t) = f2{Xs,P,t) 



P 
Pin] 



■ exp [-qC{z)] 



(13) 



where /l^'^'^(p) is given in Equation ([7j) and z = p/p*. The function C{z) is defined as 

dz' 1 



z' exp (1/2;') 



(14) 



where Zjnj = Pinj/p* (IKang fc Ryull2010l ). Of course, for p > p*, the acceleration is limited 



by the shock age and so pre-existing CRs will be simply advected downstream, resulting in 
fs{p) ~ fo{p)- These particles, however, do not make any significant contribution to the 
downstream CR pressure, if the pre-existing power-law spectrum has the slope s > 4 (see 
below) . 
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4. COMPARISON WITH NUMERICAL SOLUTIONS 

4.1. Set-up for DSA Simulations 

We carried out kinetic DSA simulations in order to test the time-dependent features of 
the test-particle solution in Equation f|T3|) . Also for shocks with typically M > a few, the evo- 
lution of CR-modified shocks should be followed by DSA simulations, because the nonlinear 



feedback of CRs becomes important (IKang fc Rydl201Cll ). We used the CRASH (Cosmic- Ray 



Acceleration SHock) code for quasi-parallel shock, in which the diffusion-convection equation 
(IH) is solved along with the gasdynamic equation modified for the effects of the CR pressure 



(IKang et all 120021 ). 



We considered shocks with a wide range of Mach number, M = 1.5 — 5, propagating into 
typical ICMs and cluster outskirts of Tq = 10'' K; the shock speed is Us = M-474 km s~^. The 
diffusion in Equation was used. In the code units, the diffusion coefficient is normalized 
with Ko = lO^K* for numerical simulations. Then, the length and time scales are given 
as lo = i^o/us and to = i^o/u^, respectively. Since the ffow structure and Pc profile evolve 
self-similarly, a specific physical value of matters only in the determination of Pmax at a 
given simulation time. For instance, Pmax/'^pC ~ 10^ is achieved by the termination time of 
t/to = 10 in our simulations. Simulations start with purely gasdynamic shocks initially at 
rest at Xs = 0, and the gas adiabatic index is 7^ = 5/3. 

As for the pre-existing CRs, we adopted fo{p) = fprdp/Pini)''^ for their spectrum. The 
amplitude, /pre, is set by the ratio of the upstream CR to gas pressure, R = Pcfi/Pg^, and 
we consider R = 0.01 — 0.1. We note that with the same value of R, the amplitude /pre 
is larger for softer pre-existing spectrum, i.e., larger s. To examine the effects of Alfvenic 
drift, in addition to the models with 6 = 0, we consider 6 = 0.42 as a fiducial value, which 
corresponds to Eb ~ O.lEg, i.e., the magnetic field energy density of ~ 10 % of the gas 
thermal energy density. Finally, we consider = 0.25 — 0.3 for the injection parameter. 



4.2. CR Proton Spectrum and CR Pressure 

Figure 3 shows the CR pressure profile and the CR distribution at the shock location, 
fs, from DSA simulations for a Mach 3 shock. In the cases with pre-existing CRs in (b) 
and (c), the steady-state solution without injection given in Equation ([7]) (dot-dashed line) 
is also shown for comparison. As CRs are accelerated to ever high energies (p ma^ oc t), the 



scale length of the CR pressure increases linearly with time, /^(pmax) oc Ust (IKang et al. 



2OO9I ). Left panels demonstrate that the CR pressure profile evolves in a self-similar fashion, 
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depending approximately only on the similarity variable, x/{ust). Right panels indicate 
that fs can be well approximated with the form in Equation f|T3l) . i.e., the acceleration of 
pre-existing and injected CRs along with an exponential cutoff at Pmax(^)- 

Comparing the cases in (a) and (b), we see that with the same injection parameter, 
the presence of pre-existing CRs results in higher downstream CR pressure, and that the 
re-accelerated pre-existing population dominates over the injected population. The presence 
of pre-existing CRs acts effectively as a higher injection rate than the thermal leakage alone, 
leading to the greatly enhanced CR acceleration efficiency. For the case with = 0.3 in 
(c), the injection rate is much higher than that of the case with = 0.25, yet the injected 
population makes a non-negligible contribution only near pinj. 

In Figure 4, we compare the spectrum of re-accelerated CRs from the steady-state 
solutions without injection (left panels) and the CR spectrum at the shock location from 
the time-dependent solutions of DSA simulations at t/to = 10 (right panels), in order to 
demonstrate the relative importance of the acceleration of the pre-existing and the injected 
populations. Different values of M and s are considered, but R = 0.05, 6 = 0.42, and 
eB = 0.25 are fixed. As noted before, with the same R, the amplitude /pre of the pre- 
existing CR spectrum is larger for larger s, so the re-acceleration of pre-existing population 
is relatively more important. The figure indicates that for most cases considered, the re- 
accelerated pre-existing population dominates over the injected population for the considered 
range of Mach number. Only for the cases with s = 4 and M > 3, the freshly injected 
population makes a noticeable contribution. 

Figure 5 shows the downstream CR pressure, Pc,2, relative to the shock ram pressure, 
Pouj., and to the downstream gas thermal pressure, Pg^2, as a function of shock Mach number 
M for different values of R, s, and 6. Again, = 0.25 in all the cases. As shown in the 
top panels, without pre-existing CRs, both Pc,2/ Poul and -Pc,2/-Pg,2 steeply increase with M, 
because both the injection and acceleration efficiencies depend strongly on M. For shocks 
with M > 5, Pc,2/{poul) > 0.1 and the nonlinear feedback begins to be noticeable. The 
feedback reduces the CR injection and saturates the CR acceleration, so Pc,2 from DSA 
simulations becom es smaller than the analytic estimates in the test-particle limit (see also 



Kang &: Ryull2010l ). Also the top panels compare the models with 5 = and 6 = 0.42, 
demonstrating that the Alfvenic drift softens the accelerated spectrum and reduces the CR 
pressure. 

In (b) panels, the cases with different upstream CR pressure fractions are compared: 
Pc,2 increases almost linearly with R at shocks with M < 3 in the test-particle regime, while 
the CR acceleration begins to show the saturation effect for M > 4. With pre-existing 
CRs, both Pc,2/ Pouj, and Pc,2/Pg,2 are substantially larger, compared to the case with R = 0, 
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especially for M < 3, confirming the dominance of the re-accelerated pre-existing population 
over the injected population at weak shocks. 

In (c) panels, the cases with different pre-existing slopes are compared; with softer 
spectrum (larger s), the amplitude /pre is larger and the CR acceleration is more efficient, as 
described above with Figure 4. In (d) panels, the same cases as in (c) panels except 6 = 
are shown, demonstrating again the effects of Alfvenic drift. 

These results indicate that at shocks with M < 3 in ICMs and cluster outskirts, the 
downstream CR pressure is typically a few to 10 % of either the shock ram pressure or the 
downstream gas thermal pressure. Even in the cases where the pre-existing CR population 
takes up to 10 % of the gas thermal pressure in the upstream flow, Pr.pJPg,'?. S. 0-1 in 



the downstream flo w. This is consistent with the Fermi upper limit (lAbdo et al.l 12010 



Donnert et al.ll2010[ ). 



5. CR ELECTRONS AND SYNCHROTRON RADIATION 

Since DSA operates on relativistic particles of the same rigidity {R = pc/Ze) in the same 
way, both electrons and protons are expected to be accelerated at shocks. However, electrons 
lose energy, mainly by synchrotron emission and Inverse Compton (IC) scattering, and the 
injection of postshock thermal electrons is believed to be much less efficient, compared to 
protons. 

The maximum energy of CR electrons accelerated at shocks can be estimated by the 
condition that the momentum gain per cycle by DSA is equal to the synchrotron /IC los s 



per cycle, i.e., (Ap)dsa + (Ap)rad=0 (see lWebb et al.lll984j : IZirakashvili fc Aharonianll2007l ). 



With the assumed Bohm-type diffusion coefficient, the electron spectrum has a cutoff at 



mlc^ Us / Bo ,. - 



C VlO^kmS-l-V V^[(^0.^ffA'"G;)2 + (i?2,efl/l/^G)2]' ^ ^ 

where B^s = {B^ + ^^mb)^^^ with Scmb = 3.24 x lO"*' G is the effective magnetic field 
strength for synchrotron and IC coolings upstream and downstream of shock, and 6 = was 
assumed. Note that the electron cutoff energy is a time-asymptotic quantity that depends 
only on the shock speed and the magnetic field strength, independent of the shock age. 
For a Mach 3 shock and Bq = 1 fiG, for example, the shock jump condition gives cr = 3, 
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q = 4.5 (with 5 = 0) and B2 = 3/iG (assuming B cc p), resulting in the cutoff Lorentz factor, 
7e,cut = Pcnt/rrieC ^ 5.6 X 10^ (Ms/1000 km s~^). 



Thus, we may model the downstream electron spectrum as 



p 



Pent 



(16) 



where fp,2{p) is the downstream proton spectrum (jZirakashvili fc AharonianI 120071 ). The 
elect ron-to-proton n umber ratio, K^/p, is not yet constrained precisely by plasma physics (see. 



e^ 



Reynolds! I2OO8I ). Although K^/p ~ 10 ^ is inferred for the Galactic CRs (ISchlickeiser 
2002h. a much sma ller value, Kf,/p < 10"'^, is preferred for young supernova remnants 



(IMorlino et al.ll2009l ). However, Ke/p for the pre-existing population in ICMs and cluster 
outskirts could be quite different from these estimates. 

Next, from the electron spectrum in Equation flTB]) . we consider the synchrotron emission. 
The averaged rate of synchrotron emission at photon frequency u from a single relativistic 
electron with Lorentz factor 7e can be written as 



(P,(7e)) = -CaTf3'UBl!Mle), 



(17) 



where /3 is the particle speed in u nits of c, a r is the Thomson cross section, and Ub is the 
magnetic energy density (see, e.g., I1991I ). The frequency distribution function, 0j^(7e), 
which satisfies the normalization J (f)y['^)dv = 1, peaks at 



-^peak 



280 



B 



1 /iG 



_7£ 

10^ 



MHz, 



:i^ 



where ul = eB/rUeC is the Larmor frequency. If we approximate that the synchrotron 
radiation is emitted mostly at z/ = z/peak {i-e.., (pu{l) is replaced by a delta function centered 
at z/ = z/peak), the synchrotron volume emissivity from the CR electron number density, 
ne{'je)d'je = fe{p)p'^dp, becomes 



J(z/) ^ ^ccrT/3^?7ij— ne(7e 
6 Vl 



(19) 



with 7e corresponding to the given z/peak = z/ in Equation f|T8|) . So the ratio of the downstream 
to upstream synchrotron emissivity at a given frequency v can be written as 



J2{t^) _ B2 7e,2/e,2(7e,2) 
Jo(z/) 5o7fn/e,o(7e,o)' 



(20) 



where 7e,o and 7e,2 are the Lorenz factor that corresponds to the given z/peak = z/ in Equation 
ffTSjl for upstream field Bq and downstream field B21 respectively. 
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For power-law spectra, the ratio J2{'^)/Jo{^) can be written in a more intuitive form. 
If the ratio Kf./p of the pre-existing population is comparable to or greater than that of 
the injected population, pre-existing electrons are more important than injected electrons at 
weak shocks of M < 3, as pointed out in the previous section. Then, the downstream electron 
spectrum /e,2 can be approximated by the distribution function in Equation ([7]) with a 
Gaussian cutoff, exp(— p^/p^^J. Again adopting fe,o{le) oc 7"^* for pre-existing CR electrons, 
the downstream spectrum is fe,2{je) oc 7"^, (unless q = s) for 7e < 7e,cut = Vcnt/f^eC- Then, 
the ratio of the downstream to upstream synchrotron emissivity at v becomes 



(r-l)/2^ 



/e,2(7e) 
Jefiile) 



.280 MHz 



(s-r)/2 



^«>(r-l)/2 p-(s-r)/2 
^ -°0,AtG 



/e,2(7e 



= 10* 



280 MHz 



(s-r)/2 



(21) 



where -Bq.mG and i?2,AtG are the upstream and downstream magnetic field strengths in units of 
yuG. In the second step, we assumed that -B2/-B0 = (P2/P0)"' = where w = 1 corresponds 
to B (X p implied by the diffusion model in Equation (jlj). 

Figure 6 shows /e,o(7e)//e,2(7e) at 7e = 10^, and (72/^1)230 = J2i^)/M^) at z/ = 280 
MHz for Bq = 1/iG and w = 1 for the cases considered in Figure 5. Here, we assume that 
Ke/p is the same for both the pre-existing and injected populations. Since the electron cutoff 
momentum is 7cut ~ 10^ for the shock parameters considered here, the choice of 7e = 10^ 
and u = 280 MHz (see Equation flTH]) ) as the representative values should be safe. As shown 
in Figure 5, for M < 3, the downstream CR proton pressure can absorb typically only a few 
to 10% of the shock ram pressure even for R = 0.05. Yet, the acceleration of CR electrons 
can result in a substantial enhancement in synchrotron radiation across the shock. Our 
estimation indicates that the enhancement factor, {J2/ Ji)2m-, can be up to several at shocks 
with M ~ 1.5, up to several 10s for M ~ 2, and up to several 100s for M ~ 3. This is partly 
due to the large enhancement of the electron population across the shock, /e,2//e,05 which is 
typically an order of magnitude smaller than the ratio {J2/ Jo)2sa- Additional enhancement 
comes from the amplification of magnetic fields across the shock, B2/ Bq. 

We note that for the compression of a uniform magnetic field, B oc p^^^, that is, w = 2/3. 
With this scaling, {J2/Jo)280 should be a bit smaller than that in Figure 6. However, it 
is also quite plausible that the downstream magnetic field is stronger than that expected 
for simple compression. It has been suggested that at shocks, especially at strong shocks, 



the downstrearn magnetic field is amplified by plasma instabilities (see, e.g., iLucek fc Bell 



2OOOI : lBellll2004l ). although the existence of such instabilities has not been fully explored for 



weak shocks. Moreover, the magnetic field can be further amplifie d by the turbulence that is 



induced through cascade of the vorticity generated behind shocks (iGiacalone fc Jokipiill2007 
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Ryu et al.ll2008l ). In such cases, the ratio {J2/ Jo)2so could be larger than that in Figure 6. 
In that sense, our estimate for the synchrotron enhancement factor may be considered as 
conservative one. We also note that with s > r in Equation fl?Il) . J2{i^) / Jo{i^) is larger at 
higher frequencies, but smaller with larger Bq. 

The above enhancement in synchrotron emission across the shock can be compared to 
the enhancement in Bremsstrahlung X-ray. The Bremsstrahlung X-ray emissivity is given 
as Jx oc p'^y/T, so the ratio of the downstream to upstream emissivity can be written as 

Jx,2 _ ^2/^2 



Jx,o V Tq 



4M2 \'/V5M2-i^'/' 



M2 + 3 / V 4 



(22) 



in the limit where the CR pressure does not modify the shock structure. The enhancement 
in Bremsstrahlung X-ray emission, Jx,2/Jx,0y is 3.6, 7.5, and 17 for M = 1.5, 2, and 3, 
respectively. These values are substantially smaller than {J2/ JQ)2m shown in Figure 6. This 
implies that shocks in ICMs and cluster outskirts may appear as radio relics, but not be 



detected in X-ray, for instance, as in the case of CIZA J2242. 8-1-5301 (Ivan Weeren et al. 



2OIOI). 

Since the synchrotron/IC cooling time scales as 

13 

7e \l pG 



t,. . f «fi I (23) 



( IWebb et al.lll984j ). behind the shock the width of the distribution of CR electrons with 7e 
becomes d ^ M2^rad(7e) oc 7~^. For instance, electrons radiating synchrotron at z/ ~ 1 GHz 
have mostly the Lorentz factor of 7e ~ 10^ in the magnetic field of -B2 ~ a few fiG. So 
the width of the synchrotron emitting region behind the shock is d ~ M2^rad(7e = 10^) ~ 
100 kpc (ti2/10^ km s~^) as long as the shock age t > trad(7e = 10^) ^ 10^ yrs. This is indee d 



of order the width of bright radio relics such as CIZA J2242. 8-1-5301 ( Ivan Weeren et al.ll2010l ). 



Moreover, from the fact that d oc 7^^, we can identify another feature in the inte- 
grated synchrotron spectrum. The volume integrated electron spectrum, -Fe,2(7e) oc /e,2(7e) • 
c? oc 7e steepens by one power of 7e above the break Lorentz factor, 7e^br ~ 2.45 x 
10^ (10^ yrs/t) (-Beff,2/1 pG)^^, where t is the shock age. Note that the break Lorentz factor 
is basically derived from the condition, t = trad in Equation ( 12^ and so independent of the 
shock speed. Hence, if /e,2(7e) oc 7"'', in observations of unresolved sources, the integrated 
synchrotron emission, S^, oc z/~", has the spectral slope a = (r — 3)/2 for u < z/br and 
a = (r — 2)/2 for z^br ^ ^ ^cut- Here, the two characteristic frequencies, z/br and Ucut, 
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correspond to the peak frequency in Equation ffTSjl for 7e,br and 7e,cut, respectively. So the 
spectral slope of the integrated spectrum just below the cutoff frequency is steeper by 0.5 
than that of the resolved spectrum. 

6. SUMMARY 

Cosmological shocks are expected to be present in the large-scale structure of the uni- 
verse. They form typically with Mach number up to 10^ and speed up to a few 1000 km 
s~^ at the present universe. Shocks in ICMs and cluster outskirts with relatively high X-ray 
luminosity, in particular, have the best chance to be detected, so they have started to be 
observed as X-ray shocks and radio relics (see Introduction for references). Those shocks are 
mostly weak with small Mach number of M < 3, because they form in the hot gas of kT > 
keV. 

In this paper, we have studied DSA at weak cosmological shocks. Since the test-particle 
solutions could provide a simple yet reasonable description for weak shocks, we first sug- 
gested analytic solutions which describe the time- dependent DSA in the test-particle regime, 
including both the pre-existing and injected CR populations. We adopted a thermal leakage 
injection model to emulate the acceleration of suprathermal particles into the CR population, 
along with a simple transport model in which Alfven waves self-excited by the CR streaming 
instability drift relative to the bulk plasma upstream of the gas subshock. 

We then performed kinetic DSA simulations and compared the analytic and numeri- 
cal solutions for wide ranges of model parameters relevant for shocks in ICMs and cluster 
outskirts: the shock Mach number M = 1.5 — 5, the slope of the pre-existing CR spectrum 
s = 4 — 5, the ratio of the upstream CR to gas pressure R = 0.01 — 0.1, the injection param- 
eter eB = 0.25 — 0.3, and the Alfvenic speed parameter 6 = — 0.42. The upstream gas was 
assumed to be fully ionized with Tq = 10'' K. 

The main results can be summarized as follows: 

1) For weak shocks with M < 3, the test-particle solutions given in Equation (fT3|) should 
provide a good approximation for the time-dependent CR spectrum at the shock location. 
We note that the test-particle slope, q, in Equation and the maximum momentum, 
Pmax(^), in Equation (E]) may include the Alfvenic drift effect. 

2) For the injection parameter considered here, es = 0.25 — 0.3, the injection fraction 
is rather low, typically ^ ~ 5 x 10^^ to 10^^ for M < 3. The pre-existing CR population 
provides more particles for DSA than the freshly injected population. Hence, the pre-existing 
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population dominates over the injected population. If there exist no CRs upstream [R = 0), 
the downstream CR pressure absorbs typically much less than ~ 1 % of the shock ram 
pressure for M < 3. With pre-existing CRs that accounts for 5 % of the gas thermal 
pressure in the upstream flow, the CR acceleration efficiency increases to a few to 10 % for 
those weak shocks. 

3) For the pre-exisiting population, the enhancement of the distribution function across 
the shock, f2{p)/fi{p), at a given momentum is substantially larger than that expected from 
the simple adiabatic compression. Hence, with amplified magnetic fields downstream, the 
re-acceleration of pre-existing CR electrons can result in a substantial synchrotron radiation 
behind the shock. We estimated that the enhancement in synchrotron radiation across the 
shock, J2{^)/Jo{^), is about a few to several for M ~ 1.5, while it could reach to 10^ — 10^ 
for M ~ 3, depending on the detail model parameters. This is substantially larger than the 
enhancement in X-ray emission. 

4) Unlike protons, relativistic electrons lose energy by synchrotron emission and IC 
scattering behind the shock, resulting in a finite width of synchrotron emitting region. In 
ICMs and cluster outskirts with fiG fields, the radio synchrotron emission at z/ ~ IGHz 
originate mostly from the relativistic electrons with 7e ~ 10^, which cool in a time scale of 
^rad ~ 10*^ yrs. So the width of the ~ 1 GHz synchrotron emitting region is d ~ U2trad ~ 
100 kpc (n^/lOOO km s~^) for a shock of age t > trad- 

Finally, although the CRASH numerical code and our thermal leakage model are de- 
veloped for quasi-parallel shocks, the main conclusions in this paper should be valid for 
quasi-perpendicular shocks as well. It is recognized that the injection may be less efficient 
and the self-excited waves are absent at perpendicular shocks. However, bo th of these prob- 



lems are alleviat ed in the presence of pre-existing CRs and turbulence ( iGiacalond 12005 



Zank et al.l 120061 ). So the diffusion approximation should be valid and the re-acceleration of 
pre-existing CRs are similar at both kinds of shocks. Then, we expect our results can be 
applied to, for instance, CIZA J2242. 8+5301, the radio relic whose magnetic field direction 
inferred from the polarization observation is perpendicular to the shock normal. 
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Fig. 1. — Surface area of shocks in ICMs and cluster outskirts, S (red solid line), and kinetic 
energy flux passed through surfaces of the shocks, Fgh (blue dashed line), as a function of 
Mach number M at z = 0. Only shocks with the preshock gas temperature of Ti > 10^ K 
are considered. 




Fig. 2. — Steady-state solution for the downstream CR spectrum, f2^^{p)^ given in Equation 
([7]) (sold lines), accelerated from an upstream CR spectrum, /o(p) oc p^"* (dashed lines). A 
shock with Mach number M = 3 is considered, so the test-particle slope is g = 4.5 (with 
5 = va/cs = 0). The CR injection is ignored and the distribution function f{j))p^ is plotted. 
Left: The case with the slope s = 5. The dotted lines show the time-dependent solution at 
the shock location, fs{p), from the corresponding DSA simulation. Right: The case with 
slope s = 4 and 4.5. 
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Fig. 3. — Time-dependent solution from DSA simulations of a Macli 3 shock. The CR 
pressure profile as a function of the similarity variable x/{ust) (right panels) and the CR 
distribution at the shock location, (left panels) at t/to = 0.5 (red dotted lines), 1 (blue 

dashed lines), 5 (green long dashed line), and 10 (black solid lines) for three cases. The top 
panels show the case with the injection parameter = 0.25 and without pre-existing CRs. 
The middle panels show the case with = 0.25 and with pre-existing CRs: the ratio of the 
upstream, pre-existing CR to gas pressure is R = Pcfi/Pgfi = 0.05 and the spectral slope of 
the pre-existing CRs is s = 4.5. The bottom panels show the case with the same pre-existing 
CRs, but with a higher injection rate, = 0.3. All cases shown have 6 = va/cs = 0.42. 
In the left panels, the CR pressure is displayed in different vertical scales for clarity. In the 
right panels, for the cases (b) and (c), the (magenta) dot-dashed lines show the steady-state 
solution of the re- accelerated CRs, fl^^^ip), in Equation ([7]), while the (cyan) dot-long dashed 
lines show the pre-existing CR spectrum, fo{p), for comparison. 
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log(f2p4): steady-state log(fsP^) at t/t^=10 
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Fig. 4. — Left: Steady-state solution of the re-accelerated CR spectrum (without injected 
population), fl^^^ip), in Equation ([7]). Right: Time- dependent solution of the CR distri- 
bution at the shock location, fs{p), at t/to = 10 from DSA simulations with the injection 
parameter = 0.25. Three different spectral slopes of pre-existing CRs s = 4, 4.5, 5 are 
considered. In all cases, R = Pc,o/Pg,o = 0.05 and S = va/cs = 0.42. Results are shown for 
shocks with Mach number M = 1.5 (red dotted lines), 2 (blue dashed lines), 3 (green long 
dashed lines), 4 (magenta dot-dashed line), 5 (black solid lines). The (cyan) dot-long dashed 
lines plot the pre-existing CRs, fo{p). 
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Fig. 5. — Ratios of the downstream CR pressure to the shock ram pressure (left panels) 
and to the downstream gas thermal pressure (right panels) for different model parameters. 
Lines show the ratios estimated from the analytic formula in Equation ( |T3|) . while symbols 
show the time-asymptotic values from the corresponding DSA simulations at t/to = 10. 
The top panes show the cases without pre-existing CRs for two different 6 = va/cs = 
and 0.42. The second panels from top show the cases with pre-existing CRs of different 
R = Pc,o/Pg,o = 0.01, 0.05 and 0.1. The spectral slope of pre-existing CRs is s = 4.5, and 
6 = 0.42 is adopted. The third panels from top show the cases with pre-existing CRs of 
different s = 4, 4.5 and s = 5. Other parameters are R = 0.05 and 6 = 0.42. The bottom 
panels show the same cases as the third panels except 6 = 0. In all cases, = 0.25 is used. 
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Fig. 6. — Ratios of the downstream to upstream CR electrons, /e,2(p)//e,o(p) at p = 10^ mgC 
(left panels) and the downstream to upstream synchrotron emissivity, J2{y)l Ji^iy) at z/ = 280 
MHz (right panels). The synchrotron emissivity was calculated by Equation fl20p for the 
upstream and downstream magnetic field, Bq = \ fiG, and B2 = (p2/po)-Bo, respectively. 
The same cases as in Figure 5 except the one without pre-existing CRs are shown here. See 
the caption for Figure 5 for different line and symbol types. 



